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^ , Abstract. We establish the existence of two weak coupling regime effective dynamics for 

^^ i an open quantum system of repeated interactions (vanishing strength and individual inter- 

action duration, respectively). This generalizes known results [3j in that the von Neumann 
algebras describing the system and the chain element may not be of finite type. Then 
(but now assuming that the small system is of finite type), we prove that both effective 
dynamics capture the long-term behavior of the system: existence of a unique asymptotic 

(-H ' State for them implies the same property for the respective exact dynamics — provided that 

f^ the perturbation parameter is sufficiently small. The zero-th order term in a power series 

I , expansion in the perturbation parameter of such an asymptotic state is given by the asymp- 

f| ■ totic state of the eft'ective dynamics. We conclude by working out the case in which the 

"*^ ' small system and the chain element are spins. 
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1. Introduction 



Recall that an open quantum system consists of a so-called small system ,5^ immersed 

in a reservoir S^,, and that one is usually interested (perhaps by necessity) only in the 

observables of 5^ . In the repeated interaction model one assumes that the reservoir is an 

infinite chain of identical subsystems {S,^n&i, called chain elements, which interact with 

(^ ■ -5^ sequentially, one at a time, in the order given by their labels n e N. Here we will 

OO ! suppose that: 

• The time that S^ spends interacting with each S„ — which could depend on n or 
even be random — is actually constant, equal to r > 0. 
is^ \ • The way in which S^ interacts with each Sn is also independent of n. 

''-^ ■ • All chain elements are initially in the same state. 

More general models can be considered, as in 18] |2l ■ 

Repeated interaction systems (RISs) have been used in connection with several do- 
mains, including quantum optics ITSJI (in particular, regarding quantum state preparation 
lfT6]l ) and quantum noises 12J|4]|5l . From an open systems point of view, they are interesting 
because of their mixture of simplicity — they have, by construction, a markovian nature — 
and thermodynamical non-triviality. Since not much is known about statistical physics far 
from equilibrium, that makes them a promising source of examples and inspiration; nev- 
ertheless, their rigorous study is just in its beginnings. In this article, we focus on their 
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perturbative analysis: we address the question of existence of van Hove effective dynamics 
and its use in studying the eventual asymptotic states, as we explain in what follows. 

To place things in context, let us recall some known results about open systems with 
time-independent hamiltonian. In general, the evolution restricted to the small system 
satisfies a complicated integro-diff'erential equation, and one is interested in finding as- 
ymptotic regimes in which the resulting effective dynamics is simpler. One possibility is to 
assume that the coupling between the small system and its environment is small, in which 
case one must rescale time so as to see the eff'ects of the interaction: the dynamics is, then, 
composed of a fast part coming from the free evolution, and a slow part coming from the 
interaction. As it turns out, those dynamics decouple in the limit: the slow part, called van 
Hove limit, becomes markovian; the fast one becomes noise, which is the reason why the 
weak coupling regime is also called stochastic limit [T]. The mathematical study of the 
van Hove limit was begun by Davies 1101 in 1974. The fact that the slow dynamics exists 
(at least in some cases) can be seen as one justification for the use of master equations 
when studying open systems. The procedure which gives the generator of the effective 
dynamics can be understood as a dynamical Fermi golden rule; see 1121 for an exposition 
of the subject. An interesting, somewhat unexpected result is the following: if the original 
system has an asymptotic state, then it is well approximated by the asymptotic state of its 
van Hove limit. Additional information on the subject can be found in iMl . 

The study of weak coupling regimes in the case of RISs was begun by Attal and Joye 
||3l. As we will see later, there are at least two such regimes in this context: calling A the 
strength of the interacion, one has the cases A ^ 0, and r — > as A^t —> Q. In [3|, the 
existence of the slow dynamics is established for both regimes, under the hypothesis that 
both the small system and the chain element are finite-dimensional. They also study a third 
regime (t — » while A^t is kept constant) which is not perturbative anymore; it has the 
interesting feature that one can always adjust the model in such a way that the effective 
dynamics is generated by any prescribed Lindbladian. 

Our objective in this article is two-fold: 

• To generalise the results in 13] to the infinite-dimensional case. 

• To study the extent to which the previously described relation between asymptotic 
states of a given system and its van Hove limit holds for RISs. 

The precise meaning of asymptotic state in this context is provided by Bruneau, Joye and 
Merkli (9] who have proved, assuming that the small system is finite-dimensional and un- 
der an ergodicity hypothesis, that any given initial state of the small system converges, 
when t ^ oo, towards a unique time-periodic state. It is to be noted that this is not a state 
of thermal equilibrium, to start with because it is not constant, but above all because it has 
a non-vanishing entropy production; this justifies the claim above about the thermodynam- 
ical non-triviality of RISs. 

2. Mathematical setup 

Let Ms and Me be two von Neumann algebras, meant to describe the small system and 
one individual chain element. Let a'^ : Ms -^ Ms and a^ : Me -^ Me be the *-weakly- 
continuous groups of automorphisms which correspond to their free evolutions. We will 
suppose that Ms and Me are mutually commuting subalgebras of a larger von Neumann 
algebra M which is generated by themU This permits to extend a^ , a^ : M ^ M; we 
denote the derivations which generate these extended groups by Ss and 6e, respectively. 



'This amounts to identifying Ms = M5 ® 1^, Me = 15® Me and letting M = Ms « Me- 
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and we denote a'^a'^- simply by a'^^. We write £5 and Se for the set of normal states of 
Ms and Me, respectively. 

Given a self-adjoint element v e M, consider the perturbed dynamics ip'^^ generated by 
the derivation 6s + 6e + i/l[v, ■]. It is explicitely given by the convergent series 

(1) <p'sE = {Mm + Yj^.^^^^'fi'sE.kWsE^ 

k>\ 

where the i^^^^ are given by the *-weakly-convergent integrals 

(2) ip'sj^j^ ^ \ dh-- \ dh a'^^[v, ■]a-^'^ ■ ■ ■ a'^^[v, ■]q'^^. 

Jo Jo 

We are interested in the repeated interaction evolution restricted to the small system, under 
the assumption that all chain elements are initially in the/?-KMS state oje ^ £>e- Therefore, 
we consider 

(3) fL = (Es flEf^s 'p'^eIms -Ms-^Ms, 

where n E N, f 1 € [0, t[, f = nr + f 1 and Es : M — > Ms is the conditional expectation given 
by 

(4) Es (xs xe) = xs coe(xe), Vx^ eMs,XE eMe- 

Remarkl.l. The existence of Es follows from the fact that, under the isomorphisms M^^ = 
Ms ® Ifi, Me = Is ® Me and M = Ms ® Me, it can be written as the composition 

Ms ® Me > Ms ® C = Ms > Ms ® Me- 

Equation (O defines a *-weakly continuous family of completely positive maps. Ob- 
serve that the semigroup property fails, since (p[.^^ gives the correct time evolution only if we 
start at times which are integer multiples of t. Note, however, that one can define in the ob- 
vious way a two-parameter family {i/5re's'"}ri</2eK satisfying (p[^^ - (pi-i[ and lores'" 'lores' = ^res% 
for all t\ < t2 < h. This is also related to the fact that our intuitively correct formula for 
time evolution can be obtained by exponentiation of a time-dependent hamiltonian — which 
would be somewhat more rigorous. In fact, one could consider the von Neumann algebra 
which describes simultaneously the small system and the entire chain, and define there a 
hamiltonian which, depending on the instant of time, makes the small system interact with 
the adequate chain element. One would obtain a piecewise constant generator whose ex- 
ponentiation, after composition with the right conditional expectation projecting onto the 
small system, coincides with ip[^^. We will omit the simple but lengthy and notationally 
involved proof of this fact, because it does not give any insight on the problems which 
concern us in this work. For more details, see 13] 19). 

To simplify the study of the weak coupling regime, we will impose a condition on the 
perturbation v e M which ensures that there are no first order effects: 

(HI) There exists a projection po £ Me, invariant under a'^, such that 

V = pqv{\ - po) + (1 - pa)vpo. 

Remark 2.2. First order effects (as can be seen from the Dyson series) do not reflect an 
influence from the environment: they come from the part of the perturbation which can be 
interpreted as modifying the free dynamics of the small system. 

One can think of po as the projection onto the first eigenspace of Se, which could be 
interpreted as an absolute vacuum state. In this case. Hypothesis (HI) is loosely saying 
that the small system and the chain element interact only through creation and annihilation 
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processes at the chain level. To see this in more detail, we refer the reader to fS), where 
an interaction which precisely falls within this description is considered. But Hypothesis 
(HI) can perfectly apply in other, different situations, where the interpretation just given is 
not adequate. We should warn, however, against one potentially tempting interpretation: 
by GNS construction we can always assume that coe(x) = {0.e, xHe), with Qg belonging 
to a Hilbert space on which Me acts. The projection IQ^) {Q.e\ cannot take the role of po 
because it does not belong to Me- 

Proposition 2.1. The linear operator 

(5) r(i,T) = Es^5£^ eB(Ms) 

is completely positive, normal and \\T{A,t)\\ — 1. Moreover, given t > 0, the map A i— > 
T(A, t) is analytic and, if the hypothesis (HI) holds, it is also even. 

Proof. The convergence of the Dyson series shows that /I e R i-> i^J^ € B(M) is analytic; 
it follows that A i-» T{A, t) is analytic too, since 

F € B{M) M EsF\ms e B(Ms) 

is linear and bounded (observe that E5 , being a conditional expectation, has norm 1). Com- 
plete positivity and normality are a consequence of the fact that Ej and ip'^^ have these 
properties. Since T{A,t)1 = 1, by general properties of completely positive maps we also 
have that ||r(^T)|| = 1. 

Let us check the parity. Under the hypothesis (HI), the invariance of po under the free 
evolution a'^ implies — thanks to the KMS condition — that 

Es(poXsXe) = (De(poXe)xs 

= a)E{xEa'^(po))xs 

= a)E(xEPo)xs 

= Es (xs xePo), "ixs eMs,XEeME. 

Hence, 

Esix) = Esipoxpo) + Es((l - po)x{\ - po)), V.x e M. 

Using this, all we have to do is prove that, for all odd k and xs e Ms , 

POflE,k(^s)PO = (1 - PQ)'PlE,k(^s)i^ -po) = 0, 

where (p'sEk i^ defined in (|2]). But this follows again from the invariance of po and the 
relations 

Po[v,xs] ^ pQVXsil - Po) - pt)Xsv(l - Po) 

= [v,xs](l - po), 
[v,xs]po = (1 - po)[v,xs], 
which are a consequence of the fact that po and xs commute. n 

3. Van Hove LIMIT 
Schematically, we are concerned with the study of an operator of the form 

(6) (pe^(A+^B)p)« ^ [(1 + TO(A\)e'''^)f, 

where f is a projection, A the generator of a group of isometrics, B a perturbation and 
« e N. Note that the parameter that determines the perturbative nature of a given regime is 
A^t; thus, we can immediately identify three different perturbative regimes: 
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(1) T is kept constant, in which case A must go to zero. 

(2) T ^ 0. Now, A can go to zero, remain bounded or even diverge — ^provided A^t — » 
0. 

(3) T — > oo and A^t —^ 0. 

In this article we treat the first two cases. The third one, which is a priori out of the reach 
of our method, seems to oscilate with r (the example of Section|5]gives some evidence of 
this). 

To identify the adequate time scale of an effective dynamics in each of these regimes, 
note that the approximation ^ is likely to become useless when n ^ 1/{A^t^) — that is, 
when t - nr ^ 1/(/1^t). Therefore, the appropriate time scale should be s - A^rt, irrespec- 
tive of the perturbative regime which is being considered. 

3.1. A preliminary result. Here we state a simple generalization of a theorem by Davies 
ifTTI . which is an abstract weak coupling dynamics existence result. 

Theorem 3.1. Let X be a Banach space, Aq : domAo <z X ^ X the generator of a strongly 
continuous group of isometrics and Ai : R — » B(X) a norm-continuous map. Suppose that 

1 r^ 

s-lim- dte'^^Ame-"^" 

T^oa T Jo 

exists and denote it by A i(Q)\ Then, defining A{s) = Aq + eAi(e), we have that 
lim sup ||(e'"*('=>^'^e-'"*°/'= -£■''*'("'* )x||=0, 

for any sq > and x e X. 

Proof. Davies proved this result when Ai(e) is actually constant; we will get the general 
case as a consequence, by showing that 

lim sup ||e"'*(''>/'' - e'^(Ao+sAm>/E]^^ ^ q 

'^^0 se[0,.5o] 

and using the triangle inequality. By Duhamel's formula, 

dsi e'"'^^'>''(Ai{E)-Aims^'''"^^^°^"^'^°^^''- 



g.sA(E)/£ _ gS(Ao+£A,(0))/e 



Jo 



Jo 
Now, apply the Dyson expansion and use Remark p\.ll to get the estimate 

||e.>lAW/£J| ^ j|g.„(A„/£+A,(£))|| < 1 + ^ ^||Ai(£)||\ 

k>i ^■ 
which by continuity is bounded uniformly in s. Similar considerations apply to 

||g(.!-Jl)(Ao+£Ai(0))/£|| 

from which the claim follows. 
Remark 3.1. The strong limit 



1 r^ 

Ai(O)'' = s-lim- £/fe''*Mi(0)e" 
r^oo T Jo 



is the so-called spectral averaging of Ai(0) with respect to the spectrum of Aq. There are 
at least two known conditions which ensure its existence 1121 . namely: 

(1) Ao admits a total set of eigenvectors, and 



RODRIGO VARGAS 



(2) Ai(0) is compact and X is a Hilbert space. 
In the first case, Ai(O)'' is equal to 



Y^PnAmPn, 



where the f ,,'s are the spectral projections of Aq and the sum converges strongly. Observe 
that it is, in a sense, the part of Ai(0) which commutes with Aq — and this interpretation 
holds whenever the strong limit Ai(O)'' exists. 

3.2. The regime A ^> Q. To use Theorem 13. II in the repeated interaction case we start by 
restricting our attention to the discrete semigroup consisting of integer powers of T{A, t); 
otherwise said, we regard only times which are integer multiples of t. The only problem 
then is to "interpolate" the semigroup {T{A, t)")„(=n to continuous time. 

Theorem 3.2. Suppose that Hypothesis (HI) holds, as well as 

(H2a) The spectrum ofa^^ is not dense in the circle 5 ' C C. 

Let r Q C be a curve with degCF, 0) = which encircles the spectrum of a'l, choose a 
branch of logarithm analytic in the interior of Y, and define 



2ni Jr 



dz log(z)(z - aj ) ' . 



Assume, finally, that 

(H3a) {EsiflEii^ ^ s-liml- f dt e'^"¥,siplE2e 
T^oo T Jo 

Then, the norm- continuous contraction semigroup 



'^° exists. 



satisfies 



^l - e-^^^^^s..J' : Ms ^ Ms 



lim sup \\iT{A,T)^''^'"'>^a-/^'''-''''^^-<p')x\\=0, 



for all So > 0. Here, [-J denotes the integer part of its argument. 

Proof. Recall that T{A, r) = E^^J^I , whence T(0, r) = aj and there exists an e > such 
that the curve r encircles the spectrum of r(/l, r), for all /l^ < s. Define A : ]-e,s[ — > B(M) 
by 



A(A^) = ^ J* logWfe - T(A,T)r\ 



which gives an analytic function since the dependence of T in A is quadratic. Assuming 
that 



1 r 

A'(O)'' = s-hm- dt e''^M'(0)e"''*° 

T^oc T Jo 



exists. Theorem [3T| would provide the conclusion with (f^^^ - e*'^ *"' . Therefore, we have 
to prove that A'(O)'' exists and is equal to -(Esi/jJ^j)''- To do that, recall that 

T{A, r) = aj - A^Estpl^^i^l + ^(A^) 
-. al + A^T2 + O(A^). 
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Hence, 

A'(0)^ = s-limi r dt:^ f dz log(z)(z- alT'e'^oTie-'^'iz- air' 
= ^ J* log(z)(z - air' {sAjm^ J dt e'^»r2e'^°)(z - aj)"' 
^T^^^JdzlogizXz-air^ 
= r^^ log(z) 

The integration order can be reversed, since the integrand and the domain of integration 
are both bounded; the same argument justifies the exchange of strong limit and complex 
integral. Note that Hypothesis (H3a) ensures the existence of the limit. n 

Remark 3.2. The spectral projections of Aq (which is always bounded) do not necessarily 
coincide with those of 6s , so that (f^^^. and a'^ do not necessarily commute. An extreme 
case of this would be a harmonic oscilator with energy spectrum {Inn/T : n e N}. Then, if 
we take log re" = log r + \t with t e ]-n, n{, we get Aq - 0. 

Remark 3.3. In |I3|, Attal and Joye prove Theorem |3.2| when the Hilbert spaces "T/j and 'He 
upon which Mg and Me act, respectively, are finite dimensional. Their method consists in 
solving explicitely the equation 

(7) T{A, r) = 6^*'*°+''''*'^ + C>(/), 

where Aq and Ai are the unknowns. Our method, although conceptually simpler, is essen- 
tially the same. Note that the use of a logarithm makes things easier but does not provide 
an optimal result, since in infinite dimension it might be possible that equation (jT) admits 
a solution, even if the spectrum of a^ is dense in the unit circle. 



Theorem 13.21 actually allows one to understand the behavior of ip[^^ for /I «; 1 and 
arbitrary t < l/A^; in other words, the restriction to times which are integer multiples of r 
is immaterial. 

Corollary 3.3. Under the same hypothesis of Theorem 13.21 the contraction semi-group 

ifi''„ : Ms — > Ms satisfies also 

lim sup \\{<plifa-/''' - <pljf)x\\ = 0, V,vo > 0. 
Proof. Indeed, writing s/A^ - m + ti with n - \_sI{A'^t)\, one has 



\m 



sl-f -s/A^ s l\rr/ 1 \mc 'l -(nT+/|) j I 

■es ffs - V'eff = \\T(At) Es<fil^a^ - <pj\ 



•^s 



< \\T(A, rra-/' - <flg\\ + \\T{A, r)''Es(^^£a," - IdM)^^ 

The first term is controlled by Theorem 13.21 while, using the Dyson expansion. Remark 
lA.ll and the fact that Es'p'se i ~ ^' ^^^ second is bounded by 



isi<p'lE-"sE')\\^f2i^r)Ah- 
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3.3. The regime t —> 0, A^t —^ 0. This regime is, analitically, somewhat more dehcate, 
because one has to control the dependence in r of the error as ??t -^ 0. That prevents us 
from just using functional calculus as in the previous subsection. In f3l, Attal and Joye 
use a refined, but finite dimensional, version of Theorem 13. II to deal with this; however, 
their proof cannot be easily extended to the infinite dimensional case. We take a different 
approach, which consists essentially in regrouping the error terms so that one can apply 
Theorem l3. ll directlv. 

Lemma 3.4. Given constants Aq, tq > 0, we say that A and r are admissible if 

Te[0,To], a\&[Q,aItqI 

Suppose that there exists some Co > such that, for all admissible A and t, 

||r(i,T)-{l +/IVAi)q'J|| <CaA^T', 

where A\ e B{Ms) is such that t i-^ a'^Aia^' is norm continuous. Then, again for all 
admissible A and t. 



sup \\ip; 

0<s<.so 



•'/('l''-) _ ^^fc+-l'TA,)/(i-r)|| 



0(t), 



where ^o > is arbitrary. 

Proof. Thanks to the Dyson series, with s - A^t in the notation of AppendixlAl 

r(^,r)-e^('*''+'''^'*>) 

= ji^T^Ai - A^T I dt a'sAiaAal + E(A, t), 

where, using the function /2 defined in (O, 

\\E(A,t)\\ < CqA^t^ + f2(A^T^)(A\fT^ = 0(A^t\ 
Moreover, by continuity of f i-> a'^A^a'^', 

dta'^A^ag' - tAi + O(t^) 

and we conclude that, for all admissible A and t, 

|k(/l,r)-e 



Jo 



rfe+.l^rA,) <c,^2^3. 



where the constant C\ depends only on Co, Aq, tq and Ai. Now, a standard telescope 
expansion shows that 



T(A, t)'" - [e 



T(6s+i^TAi)yn 



< m e 



Y^T(A,t)''-\T(A,t) - e''^^ 



iSs+A^TA,)\r^T(Ss+A^TA,n'"-k 



■rids+A-TAi)\\ Q(X^-j3) 



But we also have, this time using the Dyson series with e - At, that 



^t(6s+A-tA,)\\ 



2_2x 



<1+C2AV, C2=/l(.lV), 



whence 



T(A, t)'" - [e 



T{Ss+A^TAi)-\m 



<C,5or(l+C2^Vr^(''^>" 
<Ci*oTe-'»'°°(i+^^'*'^'>/('*^''=0(T). 
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We conclude by observing that, writing sjiA^T) - im + t\ with m - \_sl{ATf-\, 

rres L^ J 

< ||r(/l,T)"'(Q-5' -e''(*^+'*'^'*'>)|| + ||(r(i,T)'"-e'"^<*^+-^'^'*'>)e''<*^^'*'^'*'| 

< fi(A^T^)A\^ + 0(t){1 + fi(A\^)A\^), 
which is of 0{t) for all admissible A and t. 

Theorem 3.5. Suppose that Hypothesis (HI) holds, as well as 

exists. 



(H2b) (Es[v,f)^ = 5-//mi r dta'sEs[v,-fas'. 

T^°= 1 Jo 



(H3b) f e R I— > a'^^[v, -^a^'^ e B(M) is norm continuous. 
Let Tn,An > be two sequences such that t„ —> Q,Af^T„ — > 0. Then, the semigroup 
ipl^ = e-i(^s[v,-f)X . j^^ ^ j^^ satisfies 

hm sup (p, ((^,,, " a^ - ^effMWx.M., = 0' 

n-joo o<.5<i„ 

IWI<1 
for all fixed p e {Ms), and sq > 0. 

Proof. Observe, first, that by continuity of f i-> ^^^[v, -Ja^^ one has 



Ps£,2 = I dt2 \ dtia^j,[v,-]a'^^'[v,-]a- 
Jo Jo 



since linear operator composition B(M) x B(M) — » B(M) is norm-continuous. Therefore, 
using Dyson's expansion and the evenness of T(A, t) in A, one finds that 

2 

TiA,T) = Es{l + (ii)'(y [V, -f + 0(r3))}a5£ + ()(//) 

= {l-^Es[v,f^+0(^V), 

where we have used the fact that Esff^^ - a'^Es and 0(A*t*) is, actually, (9(/1^t-') when 
A < Ct '^^. To apply Lemma [34l we have to check that 

t ^ a'sEslv, -fa-' = Esia'^^lv, ■]Q's£)^Es 

is continuous, which is direct by hypothesis. 

To conclude we would like to use Theorem |3.1| but the group 

is only *-weakly-continuous; we have to show that it admits a predual, which then by 
definition would be strongly continuous. But we know that 6s admits a predual (the 
generator of the strongly continuous group (a^)*), and therefore it suffices to see that 
(E5[v, -J^)* : M^ — » M^ leaves the sub-space of ultraweakly continuous forms invariant. 
Now, for that it is enough that E^ [v, -J^ : Ms — > Ms be ultraweak-ultraweak continuous, 
and, since E^ is positive and normal, all we have to do is prove that the operations Ms — > M 
of left and right multiplication by elements of M are ultraweak-ultraweak continuous — 
which is an elementary property of the ultraweak topology, concluding the proof. n 
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Remark 3.4. When the small system and the chain element are finite-dimensional, the 
hypothesis on the continuity of ^^^[v', Ic^'g always holds; hence, this theorem is a gener- 
alization of the one in [13J . 

4. Asymptotic state 

In this section we will suppose that the von Neumann algebra Ms is of finite type I„ — 
that is, isomorphic to M„(C). Recall that in this case all semigroups are automatically 
norm-continuous. 

The expression "asymptotic state" in the context of quantum dynamics presupposes that 
the system is being studied in the Schrodinger picture; if we actually have a completely 
positive semigroup (f' : Ms — > Ms , the evolution of states is given by 

coi(x) - coo(ifi'(x)), xe Ms, cod eSs- 

Now, the convergence a>t > cOca for every state a>o implies the weak convergence of 

(p'{x) towards a limit P{x) which defines a linear function P : Ms — > Ms- Note that 
P{x) must be a multiple of the identity, because otherwise a)u(P(x)) would depend on coq; 
therefore. 

If' ^P, P(x) = WooWl. 

Conversely, the convergence of ip' to a rank-one projection P : Ms — > Ms (whose range 
must be CI c Ms since if'{l) = 1) implies the existence of a unique asymptotic state. 

In the case of repeated interaction systems, one must take into account the fact that the 
asymptotic state, if it exists, is, in general, r-periodic ||9l; an obvious necessary condition 
for its existence is, then, that T{A, t)" —> P(A). In the next subsection we study this situation 
from an abstract viewpoint. 

4. 1 . On the analytic perturbation tlieory of matrices. In this subsection we will sup- 
pose that T : ]-£(),£(>[ —> M„(C) is an analytic function such that 1 e specr(e) and 
||r(e)|| = 1. The classical reference for this material is 1 13]. We start with a lemma which 
lies at the heart of the section. 

Lemma 4.1. For each s e ]-eo, soL let P(e) be the spectral projection of I & spec T(e). 
Suppose that T(e)" > P(e) when < s < so- Then, 

(1) e spec P(0)T'(0)P(0). Let Q be its spectral projection. 

(2) ^(0)2 = QPiO), so PiO)Q is a projection, too. 

(3) P(0'*') — limg^o P(b) exists and is a sub-projection ofP(0)Q. 

Proof. Let T be an analytic extension of T to a complex neighbourhood of zero Q c C. 
We want to prove, in the terminology of 1 13 1, that is not a branch point of 1 e spec f (0). 
Since exceptional points are isolated, in any case we can suppose that there exist m analytic 
functions P, : Q \ ]-oo, 0] —> M„(C), which are all spectral projections of f, such that 

in 

P(0) = Yj ^'<0)- 

(=1 

Now, one of these spectral projections, say Pi -: P, must correspond to the eigenvalue 
1 e spec T{s) and must therefore be an analytic extension of P. Suppose, by contradiction, 
that is a branch point of 1 of order p - I > I. We know (see lfT3l Theorem L9]) that. 
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in this case, P admits a Laurent expansion in powers of z''''' which necessarily contains 
negative powers. However, by continuity of the norm, we have 

\\P{s)\\ = Um \\T{eT\\ < liminf ||r(fi)ir = 1. 

This means that, if we approach through the real positive axis, lim£^()+ ||P(e)|| - 1. This is 
a contradiction, and we conclude that is not a branch point of 1 . In particular, P can be 
further extended to an analytic continuation of P|]o,e„[ defined on a complex neighbourhood 
of and P(0+) exists. 

Making use of P{z), each ^o 6 P(0^)C" yields an analytic choice ^(z) = P{z)^o of 
eigenvectors of T(z) with eigenvalue 1. Now, the first order term in z in the equation 
f(z)m = m is 

r(0)f (0) + r'(0)^o = r(0) 

which, pre-multiplied by P(0), gives P(0)T'(0)P(0)^q = 0. In particular, 

E spec P{0)T'(0)P{0). 

Let Q be its spectral projection. This means that Q^o = ^o for all ^o e P(0^)C", and 
therefore that P(0+) = 2^(0+). 

Next, we show that P(O^) = P{0^)Q. This follows from applying the same reasoning 
above to the (real) analytic function T{s)*. In fact: it satisfies the hypothesis of the lemma; 
the spectral projection of 1 e spec T(s)* is Pie)*; and we have that 



d 
ds 



Tie)* = T'iOy. 

£=0 



Therefore, we conclude that P(0+)* = g*P(0+)*. 

Finally, since Q is obtained by spectral calculus from P{Q)T'(Q)P{Q) and 

[P(0),P(0)r'(0)P(0)]=0, 

we have [Q, P{Q)] - 0. To conclude, it only remains to show that 

P{0)P{0^) = P(0+)P(0) = P(0+), 

for in that case P(0+) = P{Q^)P{Q)Q = PiQ)QP{Q^). But again, the equation P(0)P(0+) = 
P(Q^) just amounts to saying that the elements in P(0^)C" are eigenvectors of r(0) with 
eigenvalue 1, and P{Q^)P{Q) = P(O^) follows from applying the same reasoning to r(e)*. 

D 

Remark 4. L Note that, by analyticity of P, one has P(s) - P(O^) + (9(e) for e > 0. 

The next result, which has some independent interest, is an application of Lemma [43] 
relating the asymptotic states of a one-parameter semigroup and its van Hove limit. 

Proposition 4.2. Let A : R —> M„(C) be an analytic function, with A(0) = 2 cikPk anti- 
hennitic (we suppose that the a^'s are pairwise different), kerA(e) ^ {0} and ||e''*''''|| = 1. 
Suppose that 

S — *00 

with Tr Q — 1. Then, there exists an sq > such that 
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Proof. We first fix some notation: write 

spec A(e) - {fl,(e) : / e {0, . . . , m}], 

with flj : R ^ C continuous for all /. Since the null space of A{e) is non-trivial, we can 
suppose that oo = 0. We have the expansion 

a,(e) = a,(0) + s^'^'Xi + Oifp-'P'), 

where Ai is an eigenvalue of 2 PkA'(Q)Pk and p, E N is the branching order of fl,(0). 
Recall that the hypothesis e'^ ^''*'^"^^' > Q is equivalent to 

spec {y PkA'(0)Pk\ \{0}<z{AeC:ReA< 0], 

with e spec }] Pi.A'(Q)Pk being a semisimple eigenvalue (in fact, simple since TrQ - 1) 
and Q its spectral projection. Hence, there exists an sq > such that, except when / = 0, 
Rea,(e) < for all < £ < sq. This ensures that e''**'^' converges to the spectral projection 
of G specA(e), which we will call P(s). Now, we can make use of Lemma 143] with 
T(e) = e'*<'^>, obtaining that P(s) = P(0+) + 0(e) for all < e < eo- 
Finally, observe that 



d 
£=0 as 



e^^'^Po = Po^ 



£=() 

^oA'(0)Po 



[l+s Cds e--^(°)A'(0)e^'*("' + 0(£2)}e^(")Po 



d 



Hence, again thanks to Lemma BTTl P(0 ) is a sub-projection of Q. But Tr 2 = 1, so that 
P{s) = 2 H- 0(e) for all < e < e,,. n 

4.2. Application to Repeated Interaction Systems. We start with the regime /i — > 0. 

Theorem 4.3. Suppose that Ms is of finite type I„ and that the effective dynamics (p"^ given 
by Theorem \3.2\ converses towards a projection P : Ms — > Ms of rank 1. Then, there exists 
a Ao > and a T-periodic family a>'^ e Ss such that 



for all < A < Aq, CO e Ss <^nd x € Ms- Moreover, 

cj\{x) = - Tr{Pa's(x)) + 0{A^\\x\\). 

Proof. After the proof of Theorem l3.2l (whose hypothesis always hold in finite dimension), 
we can write T(A,t) - e'"'*^'*'' with A : ]-eo,eo[ -^ B{Ms) = M„2(C) analytic. Now, a 
direct application of Proposition |4.2| (recall from Proposition |2. ll that ||r(/l,T)|| - 1) gives 



T{A,t)'' > P + O(A^) -. P(A^). 

k—>oo 

Since 1 e Mj is a fixed point for T(A,t), the image of P(A^) is CI Q Ms. The result 
follows with 

aj'jx)^^TT(P(A^)<p',E(x)). 

Observe that cl)\(x) is r-periodic, for 

P(A^)cpIe(x) = P(A^)EsiplE(^sx) = P(A^)x. 

a 

Remark 4.2. The state x >-^ ^ Tr(Pa'^(x)) is also r-periodic since P commutes with cj and 
PMs^ CI. 
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Now, we State the result for the regime t — > 0, A^t — > 0. We face two extra difficulties: 
(1) T(A, t) cannot be seen as a function of e = A^t; hence, in order to use Proposition 



14.21 one has to parametrize analytically the convergences r -^ 0, A^t — > 0. 
(2) Once we acknowledge the necessity of the previous step, it still has to be shown 
that one can write T(A(e), T(e)) = e'"'''*'**''', with A analytic. 

Lemma 4.4. The function T : (A,t) e R^ h^ Ese^^^'^''^^"'''^^ e B(Ms) = M„i(0 can be 
written, for A and t small enough, as 

T(A,t) = e^'^^^'^'^K 

where F : K? —* B(Ms) is analytic in a neighbourhood of (Q, 0) and F{Q, Q) — 6s. 

Proof. Indeed, we have the convergent power series expansion 

oo n 



11=0 k={) \a\=n-k 



where the multiindex a belongs to N and \a\ - Yji=o o'i- After composing with Es, the 
terms with odd k vanish and we get 

oo [n/2J ^ 
/i-l m=0 

C„,„ = E5 ^ 5^° [V, -W,^--- 6f"-' [V, ■]6f- . 

\a\=n—2m 

Now, if A and t are small enough, the logarithm series 

N=\ " n=I m=0 

converges and gives the existence of an F which is analytic and satisfies e^^*'*''"''"' = T{A, t). 
Observe that Ci,o ^ 6s, so that F(0, 0)^6s. □ 

Theorem 4.5. Let A(e) and T(e), with e e R, be two meromorphic parametrizations of A 
and T such that 

A(efT(E) = e, T(e) > 0. 

Suppose that Ms is of finite type I„ and that the effective dynamics (p"^ given by Theorem 
\3.5\ converses, as s ^ 00, towards a projection P e B(Ms) of rank 1. Then, there exists a 
T(E)-periodic family to'^ € £s <^nd an sq > such that 

for all £ e ]— so, ^oD oj ^ &s ^^d x e Ms. Moreover, 

4W = -rr(P4(x)) + 0(e2|W|). 
n 

Proof. Let F{A^t, t) be the analytic function given by Lemma l4!4l and consider the family 
of one-parameter groups 

t ^ e.'^^'^K A{e) = F{A(efr{E), r(e)) = F{e, t{e)). 

Observe that A{e) is analytic. In order to relate e''^^'^' and ip''^^, fix i > and let m(e) - 
Li/(/l(e)T(e))2j, so that 

sl{A{sfT{e)) = m(e)T(E) + fi(e), < fi(e) < t(s). 
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From now on, we will drop the dependence in e of m, A, t and f i (this should cause no 
confusion). Write 



As e ^ 0, the first term vanishes by Theorem [33] and the fact that we are dealing with 
(finite) matrices. The second equals 

Since T{A, tTu-"" = QmTA{i:)^-mTA{(i) ^jjjj 

|LiA(e)/£g-.sA(0)/£ _ gmTA(£)g-mTA(0)|| _ |LmTA(e)/gfiA(E)g-/iA(e) _ jXg-mTA(0)|| 

< e''^"P=^=o ll'**"^!! |L'iA(£)g-f,A(£) _ -^11 ^ Q 

what we get is that 

</j^ff = lime-"**'=>^'^e--"*<°'/^ 

By uniqueness of both limits and generators of semigroups, we see that -^(Es[v, ■]^)^ - 
A'(O)''. Hence, we can apply Proposition |42] to conclude that there exists an eo > such 
that 

T{A{s), T{E)f > P + 0{s^) -. P{e^), 

k—>cc 

for all e e ] -sq, so[. The proof ends in the same way as that of Theorem |4.3l n 

Remark 4.3. Since we ask from t{s) to be analytic around 0, we can as well just assume 
that T(e) = s", with n > 1. Now, the restriction A{s)^t(s) = e on the parametrizations of A 
and T — which seems to be essential in our approach — implies that 

A(s) = e('-«>/2 

(further restricting « to be odd), showing that our theorem cannot say anything of a regime 
in which both A and r go to zero. 

Remark 4.4. Let s e ]-eo, eo[- The convergence 

shows that wg depends on the values A{s) and t(s), but does not depend on the choice of 
parametrizations. 



This last remark suggests that Theorem l4. 5 1 would be better stated without any reference 
to the parametrizations. To this effect, we could consider the set 



U {(A(s),T(s)):se]-eQ,so[}, 



admissible 
parametrizations 



where sq > depends on the parametrization. However, we lack any description ot this set 
which does not actually mention the parametrizations; this is the reason why we prefer to 
state Theoi'em|4.5las we did. 
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5. A CONCRETE EXAMPLE 

In the simplest instance of a repeated interaction system, both the small system and 
the chain element are spins. This case falls under the hypothesis of (|3], in which the 
effective dynamics for the regime /I — > is explicitely calculated (for some specific choice 
of the interaction). Also, in |9 1, explicit conditions for the existence of an asymptotic time- 
periodic state are found, and the asymptotic state itself is computed at zero-th order in A^. 
Here, we illustrate how this last result can be recovered as an application of Theorem l4.3l 

Let us specify the model. We choose the representation 

Ms ^Me^ M2(C), 'Hs^'He^ C\ 



and suppose that the free evolution of observables is given by the hamiltonians 

/..=(° °)eM„ /z,=(J OjeM,. 

As for the interaction, we take 

/O 1\ la b\ (0 0\ (a c 



ol-\c df\l oj-\b rf;^^^®^^- 

Finally, we assume that the chain is initially in thermal equilibrium at inverse temperature 
y8; that is, 

/■^oo ■^oi\ _ xm + xne^'^'^ 

To make calculations, let {eo, ei} be the canonical basis of C^ and consider the basis of 
M2(C) given by uu - |et) (e/|, with k, I, e {0, 1). We find that 

a's («oo) = Moo, a'g (mqi ) = e"'^ mqi , 

a's(uio) = e""'^Mio, a's(un) = mh, 

so that assuming that 5^0 and that e""'^ ^ e '^■^ , the spectral averaging in the formula for 
the generator of the eff'ective dynamics tp'^^ must be taken with respect to the projections 

Pq^Poo + Pu, P+^Poi, P-^Pio, 

where P^i - Tr(M^,(-))M/t/. Observe that, if t is small enough, e""'^ ^ e"'^'^ . 

Since we are interested in the asymptotic state of the effective dynamics when A ^ Q, 
we must compute the spectral projection of the kernel of 

(Jeff := -(Es^Jgj)" 

= -PoEsflEaPo - P-^svIe^P- - P+Es<pIe.2P+- 
Now, if (MoilcJeff |moi) and (Miol<5eff Nio) do not vanish, that spectral projection is, essen- 
tially, the one of -Po^efflPoM^cc)- Identifying f o^2(C) = C^ through the basis {moo, mi i }, this 
operator is the 2x2 matrix 

/(mOoI (5eff |mOo) <M00|(5eff|Mll)\ 
\(Mlll4ff|M00) <Mll|(5eff|Mll)/' 

But = 5eff(l) = 5eff(Moo + Mil ), SO that this matrix has the form ( ^^ ^''° ), with 

^0 = (mooI (Jeff |moo) , ^1 = (Mlll^efflMll). 

The spectral projection of its kernel is 

do + 6i \oi do 
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and we find that 
Si = 



-2 



1 + e-/^^ 
-2 



'\b\' 



-/j£|^|2 l -cost(E-S) .2 l-cosr(£ + 5) 
p i -cost(£-5) _^£ 21-cost(£' + 5) 



(£-5)2 



(£ + 5)2 



We are in a position to compute the asymptotic state of the weak limit. As Theorem l4.3l 
ensures, it coincides at order zero with the one of the restricted dynamics, computed in ||9J. 
As sufficient conditions for its existence we recover also the result in |(9l. 

Proposition 5.1. Suppose that S ^ and \b\~ + |c|" ^ 0, and let (p"^ be the effective 
dynamics given by Theorem \3.2\ There exists some tq > such that 



^o + ^i^'^llo ^o"" 



c^ivUx)) 



for all o) e &$> x e Ms and r < tq. 

Proof. Let x = ( ;^j° 'l°\ ) € Ms . The computations above show that — ^provided there is 
convergence — 

(^ijcoo + 5oJCii /l 0\ 



^sSc„ 



(X) 



So+6i 



1 



It remains to see, for every small enough r, that there is indeed convergence. 
With respect to the basis {mqo, mh, moi> mio}^ 



/' (^0 -<Jo ■ 

-6i 5i 

<M0ll^eff|M01> 

<Ml0l<5eff|Ml0). 



3eff 



The eigenvalues of this matrix are 

0, 6o + 6u <Moil<5eff |moi) and (Miol^eff |mio>, 

and we have to check that, except for 0, their real part is strictly negative. Since \b\^ + |c|2 ^ 
0, one has that 6q + 6i = Re(5o + 5i) < 0. As for the others, up to order t^ we have that 



Re<Moil4ff|Moi> = 



'a c\ l—c a — d 
2{\+&-PE)\\b d)\0 c 



+e 



-PE 



a b\ l—b a — d 



c d \ 



.-PE 



—b a - d\la b 
b ]\c d 

-c a — d\ a c 







b d 



2(1 + e f"^) 

+ e-^^(|a|2 + |fe|2 + |c|2 + |^|2 - 2ad)] 

< —i\b\^ + \cf-) < 0, 
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whereas 

Re(«iol4ffl«io)=2(i+e-/'^)U^ ^/\^-fl -cj \d - a -bj\c d 

\c d]\d - a -b] \d - a -c]\b dj , 



-T 



2 



,|2 , 1. |2 , . |2 , I .|2 



2(1 + Q~P'^) 

+ e-^^(|a|2 + \b\^ + |c|2 + \d\^ - 2ad)} 

_T-2 

< —i\b? + \c\') < 0. 



Appendix A. The Dyson series 

In this appendix we collect the results we need on the perturbation series known as the 
Dyson series. Proofs can be found in |6|, for example. 

Theorem A.l. Let X be a Banach space with predual X, and 

A() : doniA^) Q X —> X 

the generator of a *-weakly-continuous semigroup {S']teBL^. Consider the perturbation 
A(e) — Aq + sAi, where A\ € B{X). We have that A{s) generates a *-weakly-continuous 
semigroup too, which we will denote by {S{e)']te&^. It satisfies 

S{b)' ^S' + y e" f dt„-- f'dtiS"AiS'--''Ai---AiS'"-'"-'AiS'-'". 
^ Jo Jo 

Here, the integrals are defined pointwise in the weak-* topology and give a convergent 
series for every s > 0. 

Remark A. 1. Given any one-parameter *-weakly-continuous semigroup, there always exist 
constants M > I and/3 > such that ||5'|| < Me'^' (see |6 Proposition 3.1.3]). Hence, the 
«-th term 

S'„ ^ I dt„-- \ dti S''AiS'--''Ai ■ ■ ■ AiS '"-'"-' AiS'-'" 
Jo Jo 

in the Dyson series satisfies 

||5;i|<^M"+ie^'||Aiir. 

Therefore, the error after adding up the first « - 1 terms is bounded by 

(etf-" 



(8) ^'fe^' 2 -^M'^' \\A,\f =: e/^' f„{st)s'Y , 

k>n 

where /„ : R+ ^ R+ is a continuous and increasing function. 
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